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Abstract 

a 

The S = 1/2, nearest-neighbor, quantum Heisenberg antiferromagnet on the 

o " 

square lattice with spatially anisotropic couplings is reconsidered, with par- 



ticular attention to the following question: at T = 0, does Neel order develop 



5_i 

at infinitesimal interchain coupling, or is there a nonzero critical coupling? 
A heuristic renormalization group argument is presented which suggests that 
previous theoretical answers to that question are incorrect or at least incom- 
plete, and that the answer is not universal but rather depends on the micro- 
scopic details of the model under consideration. Numerical investigations of 
the nearest-neighbor model are carried out via zero-temperature series expan- 
sions about Ising and dimer Hamiltonians. The results are entirely consistent 
with a vanishing critical interchain coupling ratio R c ; if R c is finite, it is 
unlikely to substantially exceed 0.02. 
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I. INTRODUCTION 



Lately there has been considerable interest in one-dimensional to two-dimensional 
crossover in Luttinger liquids. The motivation is largely the interest in possible non-Fermi 
liquid behavior in two dimensional models for high-T c superconductors Jl|. While much of 
this work has focused on the Hubbard or t-J model, a simpler case to understand is that of 
the spatially anisotropic square-lattice spin-1/2 Heisenberg antiferro magnet, 

H= 2J JxSj-Sj + ^2 ^Sj-Sj, (1-1) 

(and see Fig. [l]) with < J y < 4. The essential question is whether Neel order sets in for 
infinitesimal 

R = Jy/Jx (1-2) 

or whether there is a finite minimum coupling ratio necessary for long range order (LRO). 
A subsidiary question is the nature of the magnetically disordered phase at small R, if the 
latter scenario holds true. While various analytical and numerical investigations of this 
question have appeared in the literature @-§J|, it appears to be far from settled. 

The purpose of this paper is twofold. In Sec. [Tl] we review previous analytic arguments 
for or against the existence of a finite minimum ratio. We then discuss a renormalization 
group framework within which previous arguments appear to be incorrect or incomplete. 
This approach indicates that the behavior is not universal and that for any particular model 
the question may only be answered by numerical investigation. Therefore in Sec. |II| we 
present such a numerical investigation of the Hamiltonian (|1.1|) based on high order series 
expansion about both the Neel and dimer phases. That analysis seems to indicate that if 
there is a critical ratio it is rather small. 



II. THEORETICAL ARGUMENTS 
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A. Spin- Wave Theory 



Spin-wave theory plays an important role in our understanding of antiferromagnetism. 
A systematic 1/S expansion (where S is the spin magnitude) is generally believed to give, 
in low orders, a fairly accurate estimate of the sublattice magnetization on various lattices. 
In particular, it correctly predicts the absence of long range order at T = for a one- 
dimensional antiferromagnet. This effect can be seen from the leading order correction to 
the sublattice magnetization. For the anisotropic 2D model of Eq. ( |1 . 1| ) , we have: 

d 2 k 1 



(St) = ± 



(2vr) 2 



i 



(2.1) 



where 



7k 



J j> COS ri>x J v cos ko. 



(2.2) 



Jx 4" Jy 

and the k integral runs over the Brillouin zone, \ki\ < it. For J y = 0, the integrand is 
independent of k y and the integral has logarithmic divergences at k x = and 7r indicating 
the absence of Neel order. For finite J y the integral is finite and for small R we obtain: 



(St) « ± 



S + -\nR 

7T 



(2.3) 



If we took as a criterion for the stability of the Neel phase that the 0(1) correction should be 
smaller than the leading 0(5") term, then we would conclude that Neel order breaks down 
at: 



R c w e 



-ttS 



(2.4) 



For S = 1/2, numerical evaluation of the integral in Eq. ( [2.1|) gives R c ~ .03367 0. 



B. Chain Mean Field Theory 

A standard method for treating dimensional crossover problems of this type is to use the 
known behavior of the antiferromagnetic susceptibility for the one-dimensional system and 
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to treat the couplings in the second dimension in mean field theory ||. (We refer to this as 
"chain" mean field theory to distinguish it from another mean field theory to be discussed 
later.) The Hamiltonian of Eq. ( |1 . 1| ) is replaced by the following mean field Hamiltonian for 
each chain: 

Hcmf = E ■ Si+i + My(S z ) E(-!)^f- ( 2 - 5 ) 

i i 

We now calculate (S z ) using this Hamiltonian and demand self-consistency. The mean field 
critical point is then determined by: 

l = 2R Xl (T N ), (2.6) 

where Xi is the zero-frequency antiferromagnetic susceptibility for the one-dimensional chain 
with J x = 1. For a half-integer spin Heisenberg antiferromagnet, Xi{T) diverges as 1/T, as 
T — > 0. This argument then predicts a finite Neel temperature: 

T N « J y . (2.7) 

No matter how small J y , the ground state is always ordered. On the other hand, for integer 
S, Xi(0) has a finite value of order 1/A, the Haldane gap. Hence this argument predicts 
a disordered ground state for integer S and R < R c with R c « A. This result is readily 
generalized to <i-dimensional systems. The factor of 2 in Eq. fl2.6|) is simply replaced by 
the number of nearest neighbor chains. Indeed, it presumably becomes exact in the limit 
d — *■ oo. It has an obvious problem for d = 2 where the Mermin- Wagner theorem tells us 
that Neel order should be impossible at any finite temperature. Nonetheless, we might be 
tempted to believe the conclusion that it does occur at T = 0. A more sensible application 
of this mean field theory for d = 2 is to the case of a staggered intra-chain coupling at T = 
(see Sec. |111 B| , below). Eq. ( |2.6| ) then becomes: 

1 = 2/?AcmfXi(^cmf)- ( 2 - 8 ) 

Here A is the ratio of coupling on alternate links, and A£. MF is its critical value within the 
chain mean field theory. Since a staggered interaction in one dimension has scaling dimension 
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1/2, we conclude that Xi(A) oc (1 — A) -2 / 3 (up to log corrections). Now 1 — A plays a role 
roughly analogous to a finite temperature; as A — > 1, the mean field theory predicts that 
Neel order sets in at 1 — A^ MF oc R 3 ^ 2 . 

C. Renormalization Group Argument 

Consider a system of quantum chains weakly coupled to each other at T = 0. They 
may be of XY or Heisenberg symmetry. Alternatively, consider a system of classical spins 
consisting of planes which are weakly coupled to each other, at finite T. In either case 
(ignoring topological terms, for the moment), we may represent the system at long length 
scales by a 3-dimensional non-linear a-model, with action or Hamiltonian: 

S = (A A?) J d 3 x[(dj) 2 + (d y $) 2 + R(dJ) 2 } (2.9) 

Here has unit length and either 2 or 3 components in the XY or Heisenberg case. As 
before, R is the ratio of inter-chain (or plane) coupling to in-chain (or plane) coupling. A 

— * 

is the ultraviolet cutoff; i.e., the field, has Fourier modes with \ki\ < A. (We cut off 
the momentum inside a cube.) The factor of A is inserted to make the three-dimensional 
coupling constant g dimensionless, as is usually done in formulating the 3D renormalization 
group (RG) equations. 

Note that the only trace of the quasi- 2- dimensionality is the anisotropy in the (d(p) 2 
terms. We can get rid of this by a rescaling of z by a factor of V^R, that is, we define a 
new length co-ordinate, z = z/y/~R. In momentum space, we define a new 3-component of 
momentum, k' 3 = \fRk^. The action now looks completely 3-dimensional except that the 
cut-off is no longer a 3D cube but a squat box of area A 2 and height y/RA. To complete 
the elimination of the anisotropy from the action, we reduce the length and width of the 
box to y/RA also, by integrating out higher momentum modes, the standard RG procedure. 
Because the inter-chain coupling is so weak, this may essentially be done using the 2D 
RG equations; in momentum space, (k' z ) 2 is so small (due to the small cut-off) that we 
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consider the planes to be essentially decoupled until we have lowered the 2D cutoff down 
to approximately y/~RA. The value of the 2D effective coupling, when the cut-off has been 
reduced to \/RA, then acts as the initial condition for further RG calculations: as we lower 
the cut-off still further we should use the isotropic 3D RG equations using g 2 (^/RA) as the 
initial value. In the limit R — > 0, the 2D coupling flows to its zero-cutoff fixed point, #2(0). 
See Fig. [| for a sketch of one possible flow diagram. 

In both XY and Heisenberg cases, the 3D system has ordered and disordered phases, 
separated by some critical coupling g c . (For g < g c the system is ordered.) Thus, whether 
or not the system orders for arbitrarily small R is determined by whether or not #2(0) < g c . 
In the case of Heisenberg symmetry with no topological term, #2(0) is naively infinite, so 
the 3D system is in the disordered phase for sufficiently small R. From another point 
of view, the 2D system develops a finite correlation length, £, as we reduce the cut-off. 
For sufficiently small R, \^RA < l/£. Then further renormalization using the 3D RG 
equations cannot eliminate this finite correlation length. In the XY case #2(0) is basically the 
renormalized dimensionless temperature. It has a finite value along the Kosterlitz-Thouless 
(KT) critical line, #2(0) < Qkt- The important question is whether or not g KT < g c . If 
it is, then an arbitrarily weak inter-plane coupling leads to order for all T < Tkt- In the 
other case, if g c < gxr, then there will be some special temperature, T c < Tkt, such that 
below this temperature an arbitrarily weak inter-plane coupling leads to LRO but above this 
temperature there is a minimum inter-plane coupling necessary for LRO. 

Numerical simulations of the classical layered XY model indicate that gxr <3 C !■ In 
this case the present approach predicts a critical coupling, g c (R), for weakly coupled planes, 
slightly above gKT- This can be estimated as the bare coupling for which the renormalized 
coupling at scale y/~RA, g2{\fRA) m g c . Using the 2D Kosterlitz renormalization group 
equations this gives: 

g c {R)-g KT <xl/{\nR)\ (2.10) 

in agreement with estimates based on plane mean field theory || and another approach [|7||| 
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somewhat closer in spirit to the present one. 

Now let us consider the half- integer spin quantum Heisenberg chain ||. The two- 
dimensional cr-model action now has an extra topological term added to it, with topological 
angle 7T. This angle itself does not renormalize, for symmetry reasons, but it has a cru- 
cial effect on the renormalization of g 2 • There is now a finite coupling critical point, #2(0), 
shown in Fig. 0. (The strong coupling phase is spontaneously dimerized; it can be reached 
with a sufficiently strong antiferromagnetic second nearest neighbor coupling.) The three- 
dimensional RG flows are expected to be the same for both integer and half-integer spin, 
although the nature of the disordered phase may depend on whether the spin is half-integer, 
odd integer or even integer, due to a Berry's phase topological term [jT0| . We assume the 



only important effect of the topological terms is to produce a finite #2(0) critical point in 
two dimensions, as shown in Fig. |2|. The situation is then similar to the classical XY case. 
Whether or not a disordered phase occurs for weak interchain coupling depends on whether 
or not g 2 (0) > g c . 

Let us contrast this approach to either of those discussed above, in Sees. [11 A| and [11 B 



First consider the calculation of the breakdown of the Neel state in lowest order spin- wave 
theory. This can be regarded as the above calculation of the renormalization of the 2D 
coupling, using only the lowest order 2D RG equations. For large 5* we begin with a very 
small bare coupling. However, the lowest order RG equations always break down for small 
enough R. Furthermore, these equations do not distinguish integer and half-integer S. In 
cr-model language, the lowest order RG is independent of the topological term. Actually the 
RG equations are independent of it to all orders in g; its effects are exponentially small in 
g. Nonetheless, for small enough R we always renormalize into the regime where these non- 
perturbative contributions are important. Blind use of the lowest order equations essentially 
leads to the conclusion that #2(0) is infinite and hence that the system is disordered at small 
enough R. 

On the other hand consider the standard mean field argument |J. Because the 2D 
susceptibility is divergent along the whole KT critical line, or in the Heisenberg case with 
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topological term, it predicts LRO in all these cases. We argue that in order to go beyond 
a mean field treatment one should consider not the 2D susceptibility but rather the renor- 
malized 2D effective coupling. Naively if the 2D susceptibility is infinite, the 2D effective 
coupling is zero. But in fact, this is not actually the case. Since the 2D effective coupling 
is finite at zero cut-off we must consider whether or not it is below g c , leading to the above 
conclusion. 

It was suggested recently by Parola et al. |§ that, for R less than a finite R c , the long 
wavelength behavior "can be interpreted in terms of decoupled one dimensional chains." 
The implausibility of this proposal can be seen by considering a finite number of chains. In 
this case, we can analyze the scaling behavior entirely in terms of the (l + l)-dimensional 
renormalization group. The inter-chain coupling, of dimension 1, is relevant. We might then 
naively expect from standard scaling arguments that for two chains there should be a gap, 
proportional to R (up to logarithms). The existence of a gap in this case has been shown 
numerically |p"T| , |T2"f . In the case of three chains, White et al. [[12] found that the gap appears 
to vanish and this has been argued to be the case for any odd number of chains. However, 
this does not necessarily mean that the system asymptotically behaves like decoupled chains. 
Indeed, it seems more likely that the low energy states would be those of a single S = 1/2 
chain. This assertion is motivated by the behavior of a single spin S chain as a function of 
S. For integer S one expects a gap while for half-integer S one expects universal gapless 
behavior which is independent of S (i.e., always in the free boson or k = 1 Wess-Zumino- 
Witten universality class). For half- integer S the density of low-lying excitations (after 
scaling out the spin-wave velocity) does not increase with increasing S. We may expect 
similar behavior for 25" coupled spin-1/2 chains: inter-chain coupling is relevant and the 
system scales away from decoupled chains to the single chain fixed point with fewer gapless 
degrees of freedom. Note that in this case, the number of low-lying degrees of freedom does 
not scale with the area of the system but only with the length. It is of course possible that, as 
the number of chains is increased, the gap to other excitations decreases and asymptotically 
approaches zero as the number of chains goes to oo. Indeed, this must happen if the oo-chain 
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system Neel orders. However, it then seems unlikely that the low-lying excitations would 
be those of decoupled chains. We rather expect that, if a disordered phase exists at small 
R, it has a genuinely 2-dimensional nature. Previous work on the isotropic 5=1/2 square 
lattice Heisenberg antiferromagnet with next-nearest neighbor interactions suggests that a 
dimer-ordered ground state is the most likely to occur in a magnetically disordered phase 
|13| . However, other possibilities with or without a gap cannot be ruled out. 

Actually, the result of our reasoning is pretty uninformative. It is simply that we cannot 
tell from RG arguments whether or not the system orders for arbitrarily weak inter-plane 
coupling, in cases where the 2D system has an infinite correlation length. We might expect 
that the existence or non-existence of a critical R is not universal. Different realizations of 
the model may have different values of g c or #2(0), since critical temperatures (or couplings) 
are in general not universal. For instance, a different phase diagram might ensue if the chains 
are coupled with a first and second nearest neighbor inter-chain coupling, both proportional 
to R. It would thus seem that the question of the existence of a disordered phase for finite 
R in a particular model must be answered by numerical work. 

III. SERIES CALCULATIONS AND ANALYSES 

We now turn to numerical studies of weakly coupled, S = 1/2 Heisenberg chains, and in 
particular the anisotropic two-dimensional lattice of Fig. [I] and Eq. ( |1 . 1|) . It is obviously not 
possible for any numerical calculation to distinguish between Jy = 0, i.e., the persistence 
of Neel order for arbitrarily small interchain couplings, and a sufficiently small but positive 
Jy. Provided that there is no strong evidence for any particular nonzero value of Jy, the 
best one can do is argue that numerical data are consistent with a vanishing Jy and offer 
reasonable upper bounds on its value. 

Our numerical studies consist of a variety of zero-temperature series expansions (i.e., 
Rayleigh-Schrodinger perturbation theory) following the cluster-expansion techniques de- 
scribed in Ref. 14. Series expansions have several advantages over finite-size calculations in 
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two-dimensional lattice quantum many-body problems. There is no need to worry about 
cluster-shape effects — which should be of particular concern in anisotropic models such 
as the one of present interest. Given comparable computing power, series calculations can 
account for much further-range correlations than would be possible in an exact diagonaliza- 
tion calculation (see the discussion in Sec. V of Ref. 14). Finally, the fact that the series 
calculations do not directly study the model of interest, but rather yields results on a one 
(or more) parameter family of models, may allow for further informative comparisons with 
approximate analytic calculations. 

Expansions were carried out about both Ising and dimer Hamiltonians. The calculations 
and analyses will be described in detail below; here let us preview the results. The results 
of both types of expansions are consistent with J c y = 0. We believe that Jy is unlikely 
to exceed 0.02 J x . This upper bound is not much smaller than the lowest-order spin-wave 
estimate (and hence by itself cannot be taken as very strong evidence that spin-wave theory 
is qualitatively wrong); however it is notably less than the value 0.1 J x suggested by Parola 
et al. [Q based on exact diagonalization of clusters with up to 32 spins. In addition, the 
Ising expansions yield estimates for the staggered magnetization and the correlation-length 
anisotropy which are in excellent agreement with spin- wave theory for J y /J x down to 0.1. 
Finally, the dimer expansions appear to be consistent with the chain mean field theory 
of Sec. |IIB| at small J y , and thus support the proposition that Jy = 0. However, that 



conclusion must be tempered by comparison of the dimer expansions for a plane of chains 
with the corresponding calculations for a pair of chains. 

A. Ising expansions 

In order to discuss the Ising expansions, we consider a generalization of the coupled-chain 
Hamiltonian ( |1 . 1| ) in which Ising anisotropy is introduced, namely 

H= E J x [S?S? + a(SfS? + SVS])] 
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+ E J y [S?S* + a(SfS* + SVSV)}. 



(3.1) 



Henceforth we choose units of energy so that J x — 1 (and J y = R). Physical quantities 
are expanded in powers of a; we have calculated the ground state energy E g , the sublattice 
magnetization M = (So), and the correlation length anisotropy (£ y /£z) 2 which is given 
explicitly by 



The energy and magnetization series were determined to order (a 2 ) 5 , while the (£y/6r) 2 
series were only calculated to order (a 2 ) 4 . Note that each value of J y requires a separate 
calculation of the series; since it is not possible to present the complete series expansions in 
a compact format they are not displayed here, but are available as supplementary material. 
Here we discuss the analysis. 

The energy series were analyzed by direct Pade approximants. The results of five different 
Pade approximants evaluated at a = 1 are shown in Fig. [5| It is clear that they extrapolate 
smoothly between the one and two dimensional limits, that is to say between J y = and 
J y = 1. This is consistent with the expansion being convergent up to a = 1 for all J y > 0, 
which is in turn consistent with Jy = 0. However, this is quite weak evidence, as we will 
discuss in connection with the dimer expansions. Perhaps the results for E g could serve best 
as a touchstone for the quality of finite-size calculations. 

To analyze the magnetization series, we first make a change of variables originally in- 
troduced by Huse, 5 = 1 — (1 — a 2 ) 1 / 2 , which removes the square-root singularity at a = 1 
expected on the grounds of spin-wave theory. The results for five Pade approximants eval- 
uated at a = 1 are shown in Fig. f|. For J y > 0.2 the approximants are well-converged and 
appear to be in good agreement with spin- wave theory. For J y < 0.2 the convergence is 
poor, but the approximants are suggestive of Neel order for all J y and with M vanishing as 




(3.2) 



where the subscript c refers to connected correlations 



(S'oSf)c — (SqSj 2 ) (Sq)(S 2 ). 



(3.3) 
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a small power of J y . Note that the chain mean field theory implies M ~ J^ 2 as J y — > 0. In 
any case, these data provide no evidence that M — > at any particular Jy > 0. 

The correlation length anisotropy series are short, however this quantity is nonsingular 
as a — > 1 and the Pade approximants evaluated at a = 1 are extremely consistent even 
down to very small J y : see Fig. [)} For J y > 0.1 there is remarkable agreement between the 
series estimates and the spin-wave calculation presented by Parola et al. For smaller J y , 
two of the three approximants indicate £ y /£ x ~^ f° r Jy between and 0.02, and the other 
approximant is ill-behaved. These data are consistent with Jy = 0, but also with a small 
critical interchain coupling. 



B. Dimer expansions 

The dimer expansions are in the variable A for properties of the Hamiltonian H = 
Hq + XHi, where for H$ we take the columnar dimer Hamiltonian 

#o= E JxSi-Sj, (3.4) 

(ij)ev 

with T> the dimer covering of the square lattice shown in Fig. |||, and Hi is the remainder of 
the couplings, so that A = 1 corresponds to the coupled-chain model of interest, namely 

Hi = E Jx Si ■ Sj + E Jy Si ■ Sj. (3-5) 

ri-T j =e x ;{ij)<£D r i -r J =e y 

The quantities for which we have obtained series expansions, to order A 7 (which involves 
evaluation of 1041 graphs), include the ground state energy E g , moments of the antiferro- 
magnetic equal-time structure factor 

M 1 =E(So-S i )7r(r,,), (3.6) 

i 

M rr = E(So-S J )7r(r 4 )r i 2 , (3.7) 

i 

(where 7r(rj)=l if lies on the same sublattice as ro, and is —1 otherwise), and the antifer- 
romagnetic susceptibility x- 
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It should be clear from the form of the Hamiltonian that all these quantities possess 
two- variable expansions, in powers of A and A J y . Hence the coefficient of A m in any of these 
series can be expresses in terms of an order- m polynomial in J y ; and calculations of the 
seventh-order series at eight values of J y allows one to determine them at any J y by means 
of polynomial interpolation. 

Before considering the details of the dimer series analysis, let us look at the notions 
underlying this approach, and what we might expect to learn from it. Properties of the Neel- 
ordered phase are inaccessible to the dimer expansions (in contrast to the Ising expansions) ; 
what they can provide are estimates of X c (J y ), the smallest value of A for a given J y at 
which the antiferromagnetic correlation length and susceptibility diverge. If one assumes 
the simplest possible phase diagrams in the J y -\ plane, sketched in Fig. [7|, then A c (J y ) < 1 
implies that the uniformly coupled system (A = 1) exhibits Neel order for that value of 
J y . Furthermore, one may compare the series estimates of X c {J y ) with approximate analytic 
calculations for the critical line. 

Let us first briefly discuss the ground state energy series. Direct Pade approximants allow 
for estimates of E g (J y ,X = 1) which are shown in Fig. |3|. The approximants are consistent 
with each other — and with the values obtained from the Ising expansions for J y as large as 
0.4. This implies that the singularity in E g at the critical line is very weak; the alternative, 
that A c > 1 for J y < 0.4, is ruled out by both the Ising expansions and the other dimer 
expansions. 

The terms of the Mi, M rr , and x series are all positive and increasing with order in 
A. Thus one may estimate X c {J y ) either by ratio analyses or inhomogeneous differential 
approximants (of which Dlog-Pade approximants are a special case) [[Uj . 

The estimates of X c (J y ) which result from consideration of Mi are poorly converged, so 
we will only discuss M„. and \. The results for these two sets of series are presented in 
Fig. |8]. Let us address what the various curves and points signify. First, six differential 
approximant estimates for \ c {J y ) are presented for several values of J y ; these particular 
approximants are chosen because they utilize all the terms in the series and they yield good 
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estimates of A c at both J y = and J y = 1. The missing approximants at any J y are either 
defective or (in one or two cases) off-scale. The approximants are reasonably consistent over 
the entire range of J y . Second, the thin solid and long-dashed lines are the estimates of 
A c based on ratio analysis of the three highest-order terms in the x an d M rr series. For a 
series J2n c n^n, the ratios c n /c n _i are plotted versus 1/n and pairwise linearly extrapolated 
to 1/n = 0; the intercepts yield estimates of 1/A C based on three consecutive terms of the 
series. The x ratio curve lies close to the differential approximants, which reflects the fact 
that the x series are extremely well-behaved. The ratio plot (corresponding to J y = 1) in 
Fig. 2 of Ref. 14 provides further evidence that the antiferromagnetic susceptibility series is 
better behaved than other dimer expansions. Both the ratio and differential approximant 
analyses suggest that \ c (J y ) < 1 for all J y > 0, and hence that there is no magnetically 
disordered phase along the line A = 1. 

Third, two mean field estimates of X c (J y ) are presented in the figure. The dotted-dashed 
line is based on "dimer mean field theory," |16[ that is, the zeroth- and first-order terms 



in x(A) are used to estimate A c . (This is equivalent to considering a single dimer subject 
to a staggered field which is determined, self-consistently, by its staggered magnetization.) 
That result is A^mf = 1/(1 + ^y)i and it has the remarkable (and accidental) feature that 
it yields the exact A c when J y — 0. The broad solid line is based on the chain mean field 



theory discussed in Sec. [II B| . The values of Xi(^) (t ne antiferromagnetic susceptibility for 
a single chain), which are the essential input into that mean field theory, were obtained by 
integrating a differential approximant to the seventh-order series presented in Ref. 14. The 
broad solid line ends at J y ~ 0.02 because the numerical estimates of xW are n °t reliable 
for A arbitrarily close to 1. Over most of the range of J y plotted, the dimer and chain mean 
field theories give quite close values of A c . The agreement is somewhat better than one would 
expect: in both cases the leading behavior of A c at large J y is 1/2 J y , but the next-order 
terms are — 3/4 and — 1/4 for the chain and dimer mean field theories, respectively, 
and at J y = 1 the difference is hardly negligible! At small J y the chain mean field theory is 
clearly in better agreement with the series extrapolations than the dimer mean field theory. 
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To briefly recapitulate, the dimer series for the two-dimensional lattice of coupled chains 
are entirely consistent with Jy = 0. However, there is reason to doubt the strength of this 
conclusion. Since A c (0) = 1 and A c (l) ~ 0.54, it would be entirely natural to conclude 
from too-short series that \ c (J y ) interpolates smoothly between these two endpoints, even 
if the correct result is that A c does not exist for < J y < J y as would be the case if 
the scenario of Fig. 0(b) were to hold. As a partial test of the reliability of the series 
estimates of X c (J y ), we have considered the problem of two coupled chains by the dimer 
series expansion method. The chains were taken to lie parallel to the bonds in the columnar 
dimer configuration. The only differences between these calculations and those preceding are 
the set of connected clusters (and their lattice constants), and that we will present results 
for M xx rather than M rr . (In fact one can use M xx to estimate X c {J y ) for the plane of chains 
as well; the differences between the estimates based on M xx and M rr are insignificant for 
the J y of interest.) 

For two chains the only critical point in the \-J y plane is (A = 1, J y — 0) Estimates 
of X c {Jy) coming from inhomogeneous differential approximants to the x an d M xx series for 
two chains are displayed in Fig. ^|, on the same scale as in Fig. |8] for ease of comparison. 
The thick solid and thin dotted-dashed curves are the chain and dimer mean field results, 
respectively, for two chains. (One obtains these from the mean field results for planes of 
chains by the substitution 2J y — > J y .) Several points are evident upon inspection. First, 
even at very small J y one can distinguish between the estimates of A c from the differential 
approximants for the plane of chains and for two chains; they seem to approach zero with 
different slopes, on these plots. This could be taken as further evidence that, for the plane 
of chains, J y is extremely small if not vanishing. 

However, for J y < 0.2, where the approximants are well-converged, one might conclude 
that there are critical points at A ~ 1, contrary to the known behavior of coupled pairs 
of chains. We believe this line segment of "pseudocritical" points reflects the existence of 
local maxima in the correlation length at A ~ 1 for fixed, small J y , where the correlation 
length at these maxima exceeds the typical cluster length in the seventh-order calculation 
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(which is roughly 10). That such pseudocritical points should exist is entirely plausible: 
numerical studies for pairs of finite chains PJTTfl find that the gap initially decreases as J y 
is increased from zero. What is somewhat disturbing is that, for J y < 0.1, the two-chain 
approximants are as consistent with chain mean field theory (for two chains) as the plane- 
of-chain approximants are with chain mean field theory (for a plane of chains). Thus it is 
conceivable that for the plane of chains, the values of A c indicated by the series at sufficiently 
small J y are pseudocritical points, rather than true critical points, as well. 

To conclude, the problem of coupled chains at small J y poses significant challenges to 
numerical studies. Although the series expansions are consistent with Jy = they do not 
rule out a small but positive Jy. If the disordered phase exists, the correlation length along 
the x direction is probably large throughout that phase, and the gap is everywhere small. 
A potentially fruitful avenue for future studies would be the consideration of models which 
include further-neighbor couplings; suitable models would have larger critical interchain 
couplings within spin-wave theory than the simplest coupled Heisenberg chain model studied 
here. 
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FIGURES 

FIG. 1. The spatially anisotropic Heisenberg model. The solid and dashed segments correspond 
to interactions J x and J y . 

FIG. 2. Renormalization group flows for the half- integer spin quantum Heisenberg models. 
Here we assume 52(0) < g c so the system orders for arbitrarily small R. 

FIG. 3. Pade approximants for the ground state energy, evaluated at a = 1. The symbols are 
approximants to Ising series, while the lines (which are connecting points with J y = 0, 0.1, 0.2, 
0.3, 0.4, 0.5, and 0.75) are approximants to dimer series. 

FIG. 4. Pade approximants to the Ising expansion for the sublattice magnetization, evaluated 
at a = 1 following the change of variables described in the text. 

FIG. 5. Pade approximants to the Ising expansion for the correlation length anisotropy, evalu- 
ated at a = 1. 



FIG. 6. The columnar dimer covering T> used for the dimer expansions; see Eq. (3.5) 



FIG. 7. Sketches of the simplest plausible phase diagrams in the first quadrant of the J y -X 
plane for the Hamiltonian ( |3.5| ) assuming that (a) Jy = 0, and (b) Jy > 0. Note that units of 
energy are chosen so that J x = 1. The phase diagrams must satisfy two contraints that follow 
from the invariance of correlation functions with respect to multiplication of the Hamiltonian by 
a constant: (J y , A) = (a, 1) and (1/a, 1) lie in the same phase, and so must (a, b) and (ab, 1/6). In 
both cases the hatched regions constitute the Neel ordered phase, which are surrounded by lines of 
critical points. The one-dimensional Heisenberg critical point at (J y , A) = (0, 1) is indicated by the 
large dot. The entire right boundaries of the plots, J y = 00, < A < 00, are also one-dimensional 
Heisenberg critical points. All other points are supposed to have only short-range correlations. 
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FIG. 8. Estimated values of the critical A as a function of J y for the coupled-chain model (|3.5|) , 
From the x series, differential approximants displayed are [2,4;- 1] (0)> [3,2;0] (□), and [3,3;-l] (<0); 
from the M rr series, [2,4;-l] (A), [3,3;-l] (v)j and [4,2;-l] (o). See the text for a discussion of the 
significance of the various curves. 

FIG. 9. The symbols correspond to estimates of the critical A for two chains, corresponding 
to the same approximants as the preceding figure (but for M xx rather than M rr ). The curves are 
discussed in the text. 
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